How atomic nuclei cluster 
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Nucleonic matter displays a quantum liquid structure, but in some cases finite nuclei behave 
like molecules composed of clusters of protons and neutrons. Clustering is a recurrent feature 
in light nuclei, from beryllium to nickel. For instance, in 12 C the Hoyle state, crucial for stellar 
nucleosynthesis, can be described as a nuclear molecule consisting of three alpha-particles. The 
mechanism of cluster formation, however, has not yet been fully understood. We show that the 
origin of clustering can be traced back to the depth of the confining nuclear potential. By employing 
the theoretical framework of energy density functionals that encompasses both cluster and quantum 
liquid-drop aspects of nuclei, it is shown that the depth of the potential determines the energy 
spacings between single-nucleon orbitals, the localization of the corresponding wave functions and, 
therefore, the degree of nucleonic density clustering. Relativistic functionals, in particular, are 
characterized by deep single-nucleon potentials. When compared to non-relativistic functionals 
that yield similar ground-state properties (binding energy, deformation, radii), they predict the 
occurrence of much more pronounced cluster structures. More generally, clustering is considered as 
a transitional phenomenon between crystalline and quantum liquid phases of fermionic systems. 
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The occurrence of molecular states in atomic nuclei 
and the formation of clusters of nucleons were already 
predicted in the 30's by von Weizsacker and Wheeler 
PQ [2]- Even though the description of nuclear dynam- 
ics became predominantly based on the concept of in- 
dependent nucleons in a mean-field potential, a renewed 
interest in clustering phenomena in the 60' s led to the de- 
velopment of dedicated theoretical methods [3 . Numer- 
ous experimental studies have revealed a wealth of data 
on clustering phenomena in light nuclei [4 J, and modern 
theoretical approaches use microscopic models that fully 
take into account single-nucleon degrees of freedom [5HZ] • 
Clustering gives rise to nuclear molecules. For instance, 
in 12 C the second + state the Hoyle state that plays a 
critical role in stellar nucleosynthesis, is predicted to dis- 
play a three-a structure |8j [9] . The binding energy of the 
a-particle, formed from two protons and two neutrons, is 
much larger than in other light nuclei. Cluster radioac- 
tivity [TO] , discovered in the 80's, is another manifesta- 
tion of clustering in atomic nuclei. Experimental signa- 
tures of clustering are usually indirect. Quasi- molecular 
resonances are probed by scattering one cluster on an- 
other, such as in the 12 C+ 12 C system [4j[TT], and cluster 
structures are also discernible in the breakup of nuclei. 
Evidence has been reported for the formation of clusters 
in ground and excited states of a number of a-conjugate 
nuclei [4 , that is nuclei with an equal even number of 
protons and neutrons, from 8 Be to 56 Ni. 

The mechanism of cluster formation in nuclei has not 
yet been fully understood. Deformation plays an impor- 
tant role because it removes the degeneracy of single- 
nucleon levels associated with spherical symmetry. At 
specific deformations the shell structure can restore de- 
generacies corresponding, for instance, to a 2:1 ratio of 
the large axis over the small axis of a quadrupole de- 
formed system4. Consequently, the restored degener- 



acy of deformed shell closures facilitates the formation 
of clusters. However, this is a rather qualitative explana- 
tion because clustering phenomena cannot generally be 
explained by accidental degeneracies. Clustering is an 
essential feature of many-nucleon dynamics that coex- 
ists with the nuclear mean-field. Therefore, although in 
most cluster models the existence of such structures is 
assumed a priori and the corresponding effective interac- 
tions are adjusted to the binding energies and scattering 
phase shifts of these configurations, a fully microscopic 
understanding of cluster formation necessitates a more 
general description that encompasses both cluster and 
quantum liquid-drop aspects in light as well as in heavier 
nuclei. The aim of this work is to address the origin of 
clustering, i.e. to find out the conditions for cluster for- 
mation in ground states of finite nuclei, starting from a 
fully microscopic description based on the framework of 
energy density functionals. 

At present the only comprehensive approach to nu- 
clear structure is based on the framework of energy den- 
sity functionals (EDFs). Nuclear EDFs enable a com- 
plete and accurate description of ground-state properties 
and collective excitations over the whole nuclide chart 
[T2HT4] . In practical implementations nuclear EDFs are 
analogous to Kohn-Sham Density Functional Theory, the 
most widely used method for electronic structure calcu- 
lations in condensed matter physics and quantum chem- 
istry. In the nuclear case the many-body dynamics is 
represented by independent nucleons moving in a local 
self-consistent mean-field potential that corresponds to 
the actual density and current distribution of a given nu- 
cleus. Both non-relativistic and relativistic realizations 
of EDFs are employed in studies of nuclear matter and fi- 
nite nuclei. A nuclear EDF is universal in the sense that, 
for a given inter-nucleon interaction, it has the same func- 
tional form for all systems. Using a small set of global 



parameters adjusted to empirical properties of homoge- 
neous nuclear matter and data on finite nuclei, a universal 
functional provides a description of the structure of nu- 
clei across the chart of nuclides and, therefore, describes 
the coexistence of cluster and quantum-liquid aspects of 
light nuclei. 

A number of recent studies based on nuclear EDFs or 
the mean-field approach have analysed cluster structures 
in a-conjugate nuclei [T5H2T] . In Fig. 1 we display the 
self-consistent ground-state densities of 20 Ne, calculated 
with two widely used functionals that are representative 
for the two classes of nuclear EDFs: the non-relativistic 
Skyrme SLy4 [22], and the relativistic functional DD- 
ME2 [23 . The equilibrium shape of 20 Ne is a prolate, axi- 
ally symmetric quadrupole ellipsoid. Although they have 
not been specifically adjusted to this mass region, both 
functionals reproduce the empirical ground-state prop- 
erties of this nucleus: the experimental binding energy 
160.6 MeV, the radius of the proton distribution 2.90 
fm [24], and the radius of the matter distribution 2.85 
fm [25], with a typical accuracy of 1%. It is remarkable 
that, although these functionals predict similar values for 
the binding energy, charge and matter radii, as well as 
the quadrupole deformation of equilibrium shape of 20 Ne, 
yet the corresponding single-nucleon densities are quali- 
tatively very different. The density calculated with SLy4 
displays a smooth behaviour characteristic for a Fermi 
liquid, with an extended surface region in which the den- 
sity very gradually decreases from the central value of 
0.16 fm -3 . The relativistic functional DD-ME2, on the 
other hand, predicts an equilibrium density that is much 
more localized. The formation of cluster structures is 
clearly visible, with density spikes as large as 0.2 fm -3 , 
and a much narrower surface region. 

Understanding the difference in the equilibrium densi- 
ties of 20 Ne calculated with SLy4 and DD-ME2 is the key 
to the mechanism of cluster formation in this mass region 
of a-conjugate deformed nuclei. The axially symmetric 
deformation of the nuclear mean-field removes the de- 
generacy of spherical single-nucleon levels, and nucleons 
pairwise (spin up and down) occupy orbitals character- 
ized by time-reversal degeneracy. For large deformations 
these levels can be labelled by a set of asymptotic Nils- 
son quantum numbers [26] and, because of the relatively 
weak Coulomb interaction in light nuclei, the localiza- 
tion of proton and neutron orbitals is similar in Z = N 
nuclei. In the specific case of 20 Ne, ten protons and ten 
neutrons occupy five deformed Nilsson levels, with the en- 
ergy spacing between these levels being proportional to 
the deformation of the single-nucleon potential. Figure 
2 shows the partial single-nucleon densities that corre- 
spond to highest occupied Nilsson orbital. Even without 
introducing a quantitative measure of localization, it is 
obvious that DD-ME2 predicts a much more localized 
density distribution. Similar results, i.e. more localized 
density distributions calculated with DD-ME2, are also 
obtained for the other four occupied orbitals. 

Localization of densities that correspond to single- 
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FIG. 1. Plots of self- consistent ground-state densities of 20 Ne, 
calculated with the nuclear energy density functionals: DD- 
ME2[23] (top), and Skyrme SLy4 [22 [30] (bottom). The 
densities (in units of fm -3 ) are plotted in the x-z plane of 
the intrinsic frame of reference that coincides with the princi- 
pal axes of the nucleus, with z chosen as the symmetry axis. 
The inserts show the corresponding three-dimensional den- 
sity plots and the density profiles along the symmetry axis, 
respectively. 



particle orbitals is a necessary precondition for the for- 
mation of clusters, and this effect can be traced back to 
the corresponding single-nucleon spectra. The compari- 
son of spectra calculated with the two functionals shows 
that the one obtained with DD-ME2 is more spread out, 
and the more pronounced energy spacings between single- 
particle levels are also reflected in the more localized wave 
functions and partial densities. Starting from degenerate 
spherical single-particle levels, the splitting of the corre- 
sponding Nilsson deformed states is proportional to the 
deformation and to the depth of the potential. Since the 
two functionals predict almost identical equilibrium de- 
formations and radii for 20 Ne, the different energy spac- 
ings in the single-nucleon spectra must reflect the differ- 
ence in the corresponding potentials. In fact, the self- 
consistent mean-field potential of DD-ME2 is consider- 
ably deeper than that of SLy4. In the centre of the nu- 
cleus the depth of the DD-ME2 single-neutron potential 
is -78.6 MeV, whereas the depth of the SLy4 potential is 
-69.5 MeV. The corresponding values of the single-proton 
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FIG. 2. Partial nucleon density distributions that correspond 
to highest occupied (2 protons spin up and down, and 2 neu- 
trons spin up and down) in 20 Ne: ^[NnzA] = 1/2+ [220], 
calculated with the nuclear energy density functionals: DD- 
ME2[23] (top), and SLy4[22jl30] (bottom). 



FIG. 3. Top: Harmonic oscillator potentials for three differ- 
ent values of the depth: 30, 45 and 60 MeV, with the same 
radius R = 3 fm. Bottom: the radial wave functions of the 
corresponding first p-state. The position of the maximum is 
determined by the oscillator length b. 



potentials are: -72.8 MeV for DD-ME2, and -64.6 MeV 
for SLy4. The effect of the potential depth on the local- 
ization of wave functions is schematically illustrated in 
Fig. 3 where, as an approximation to nuclear potentials, 
we plot three harmonic oscillator potentials with differ- 
ent values of the depth: 30, 45 and 60 MeV, but with the 
same radius R = 3 fm. The radial wave functions of the 
corresponding p-states are shown in the lower panel. The 
oscillator length b determines the position of the maxi- 
mum and the dispersion of the wave function [27 . The 
deeper the potential, the smaller the oscillator length, 
and the wave functions become more localised. At the 
origin of clustering is, therefore, the depth of the self- 
consistent single-nucleon mean-field potential associated 
with a given energy density functional. By performing 
a series of self-consistent mean-field calculations using 
a variety of non-relativistic and relativist ic functionals, 
not only for 20 Ne, but also for 24 Mg, 28 Si and 32 S, we 
have verified that pronounced cluster structures in de- 
formed equilibrium shapes indeed occur only for deep 
single-nucleon potentials. 

The difference of the potential depths calculated with 
DD-ME2 and SLy4 is characteristic for relativistic vs 



non-relativistic self-consistent potentials. The depth of 
a relativistic potential is determined by the difference of 
two large fields: an attractive (negative) Lorentz scalar 
potential of magnitude 400 MeV, and a repulsive Lorentz 
vector potential 320 MeV (plus the repulsive Coulomb 
potential for protons) [T2HT4] . The choice of these poten- 
tials is further constrained by the fact that their sum (~ 
700 MeV) determines the effective single-nucleon spin- 
orbit potential. In a non-relativistic approach the spin- 
orbit potential is included in a purely phenomenologi- 
cal way, with the strength of the interaction adjusted 
to empirical energy spacings between spin-orbit partner 
states. Since the relativistic scalar and vector fields de- 
termine both the effective spin-orbit potential and the 
self-consistent single-nucleon mean-field, for all relativis- 
tic functionals the latter is found to be deeper than the 
non-relativistic mean-field potentials, for which no such 
constraint arise. 

More generally, fermionic systems can exhibit a crys- 
talline phase or, on the other extreme, a quantum liquid 
phase. B. Mottelson considered the quant ality [28] pa- 
rameter to show that nuclear matter displays a quantum 
liquid structure. This concept can be generalised by con- 
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FIG. 4. Transition from crystalline to a quantum liquid phase, 
including the cluster phase, a = b/r , where b is the dis- 
persion of the fermion wave- function, and ro the typical inter- 
fermion distance. 



sidering nuclear clusters as transitional states between 
crystalline and quantum liquid phases (Fig. 4). The di- 
mensionless ratio a = b/ro , where b is the dispersion 
of the nucleon wave- function, and ro the typical inter- 
nucleon distance (~ 1.2 fm), is the natural parameter to 
quantify nuclear clustering, in analogy with similar con- 
siderations in condensed matter [29]. 
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where Vo is the depth of the potential, and R the radius 
of the system. 

When a >1 nucleons are delocalised and the nucleus 
has a quantum liquid structure. The transition to a clus- 
ter state occurs when a ~ 1, nucleons become more lo- 
calised and form a molecular structure (Fig. 4). In the 
present analysis we find that a < 1 for the relativistic 
functional, whereas a >1 for the non-relativistic func- 
tional. Moreover, from its definition in the case of a 
harmonic oscillator potential (see Eq. (1)), a obviously 
increases with the number of nucleons (nuclear radius). 
Cluster states, therefore, are less likely to appear in heav- 
ier nuclei. 
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